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Abstract. Wc consider one-dimensional asymmetric exclusion processes with a 
simple attractive interaction, where the distance between consecutive particles is not 
allowed to exceed a certain limit and investigate the consequences of this coupling on 
the transport properties in the presence of random-force type disorder by means of 
a phenomcnological random trap picture. In the phase-separated steady state of the 
model defined on a finite ring, the properties of the density profile are studied and 
the exponent governing the decay of the current with the system size in the biased 
phase is derived. In case all consecutive particles are coupled with each other and form 
a closed string, the current is found to be enhanced compared to the model without 
coupling, while if groups of consecutive particles form finite strings, the current is 
reduced. The motion of a semi-infinite string entering an initially empty lattice is 
also studied. Here, the diffusion of the head of the string is found to be anomalous, 
and two phases can be distinguished, which are characterised by different functional 
dependences of the diffusion exponent on the bias. The obtained results are checked 
by numerical simulation. 
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1. Introduction 

In low dimensions, the characteristics of transport in inhomogeneous media may 
drastically differ from those in homogeneous environment. This is to be seen already 
for the thoroughly-studied problem of random walk on a one-dimensional lattice with 
quenched random hop rates [HEIEIIIIEIEICFJIE]. The most striking anomalies are 
observed when the disorder is of random-force type, i.e. the direction of the local bias is 
random. In this case, the diffusion is anomalous, characterised by a diffusion exponent 
continuously varying with the global bias [5j [6] , while it becomes logarithmically slow for 
zero global bias — a phenomenon known as Sinai diffusion [Jj. Much less is known for the 
transport of interacting many-particle systems on disordered one-dimensional lattices. 
The zero-range process (ZRP) [9l [10], where lattice sites are allowed to be multiply 
occupied by identical particles, has a product-measure steady state even in the case 
of random hop rates and it is closely related to the corresponding one-particle problem 
(random walk) on the same lattice. Beyond ZRP, perhaps the simplest interacting driven 
many-particle system is the asymmetric simple exclusion process (ASEP)[TT|[T2]. where 
particles interact by hard-core exclusion. For this model, there are two ways to introduce 
disorder: random rates can be associated either with particles or with sites (links). 
The former model, the ASEP with particle-wise disorder [T3l [HI [15] can be mapped 
to a disordered ZRP, whereas for the ASEP with site- wise disorder pSJ [T71 [18], EHl ED], 
general exact solutions are not at our disposal. By means of a phenomenological random 
trap picture, the current of the latter model has been shown to display an anomalous 
behaviour similar to that of the random walk in the case of random-force type disorder 
PI. 

Hard-core exclusion can be thought of as a simple repulsive interaction between 
neighbouring particles. One may also consider attractive interactions between particles, 
the simplest form of which in one-dimension is, analogous to hard-core exclusion, when 
the distance between consecutive particles is not allowed to exceed a certain value, say, 
I + 1 lattice spacings (I > 1), i.e. a string of coupled particles is formed. From the point 
of view of holes (empty sites), this prescription transforms to the constraint that the 
size of clusters of holes is at most I. Regarding the number of holes in front of particle 
i as an occupation number on the ith site of a virtual lattice, the process maps to the 
so-called generalised exclusion process [11], which is a ZRP with an upper limit I for the 
local occupation numbers. Considering particle-wise disorder in the model of coupled 
particles, it transforms to a generalised exclusion process with site-wise disorder. For 
the latter model with random-force type disorder, the scaling form of the current was 
argued to be independent of /, thus identical to that of the ASEP, which is the I = 1 
limit of the generalised exclusion process [19]. The model of coupled particles with 
site-wise disorder with I = 1 maps to an exclusion process with double-sized particles 
with site- wise disorder, however, a similar mapping does not exist for / > 1 [19] . For 
the process with I — 1 in the presence of random-force type disorder, the dynamical 
exponent has been found to differ from that of the ASEP [T9] . 
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The aim of this paper is to study the exclusion process with coupled particles 
in the presence of site-wise random-force type disorder for general I and to explore, in 
which extent the transport properties are modified compared to the ASEP by the simple 
attractive interaction introduced above, which breaks particle-hole symmetry. We shall 
investigate the distribution of the sample-dependent current and the properties of the 
density profile in the non-equilibrium steady state on a finite ring, as well as non- 
stationary phenomena such as the diffusion of a semi-infinite string by means of a 
phenomenological random trap model and numerical simulation. 

Beside theoretical interest, another source of motivations is the applicability of 
simple variants of exclusion processes to the description of transport phenomena in a 
large variety of real systems such as vehicular traffic [21 J or various biological transport 
processes [22]. Concerning the latter, random-force type disorder may emerge in several 
contexts such as DNA unzipping [23], translocation of RNA or DNA through pores 
[21] or motion of molecular motors [25J on heterogeneous tracks [26, 24J. In many 
cases, molecular motors act in groups and the collective effects have been studied in 
different models [27], including models with hard core repulsion [28J or elastic coupling 
[29] between particles moving in a periodic potential. Our simple model may serve as a 
ground for testing the joint effect of hard-core repulsion, coupling and disorder in such 
model systems. A more direct example is served by a recent work, in which the diffusion 
of finite strings of coupled particles on a homogeneous one-dimensional lattice has been 
studied [30J, motivated by the modelling of synthetic molecular systems, known as 
molecular spiders, which can move on surfaces and tracks [31]. The model to be studied 
here may thus be relevant for investigating the influence of track heterogeneity on the 
transport properties of these systems. 

The rest of the paper is organised as follows. In Sec. (2], the model to be studied is 
defined. In Sec. [31 the properties of closed strings are investigated on finite rings and 
the results are formulated generally in Sec. HI In Sec. [5j the problem of diffusion of 
a semi-infinite string is discussed, while Sec. [6] is devoted to problems concerning the 
traffic of finite strings, including the steady state and the invasion. Finally, the paper 
is closed with a discussion of the results in Sec. [7J 

2. The model 

The process that we focus on in the first part of this work is defined on a ring with L sites. 
On this lattice, N particles, which are numbered consecutively along the chain from left 
to right, are distributed in such a way that each lattice site is occupied by at most one 
particle and, in case of a closed string, for the number of empty sites hi in front of the 
zth particle, hi < I holds for all i, while in case of an open string hi < I is prescribed for 
all but one particle (the Nth. one). Obviously, for the closed string, L/(l + l) < N must 
hold. In this system, a continuous-time Markov process is considered, in the course of 
which particles attempt to change their positions independently and randomly. The 
allowed transitions are the following: The zth particle on site j attempts to hop to the 
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adjacent lattice site on its right-hand side with a site-dependent rate pj, and the trial 
is successful if hi > 1 and < /. Here, h Q = and for particle 1 of an open string, 
the condition < I is cancelled. Particle i on site j hops to the adjacent site on its 
left-hand side with a site-dependent rate qj, provided > 1 and hi < I. For the iVth 
particle of an open string, the condition h^ < I is ignored again. The hop rates pj and 
qj are independent, identically distributed quenched random variables drawn from the 
distributions p{p) and 7r(g), respectively. Introducing the potential difference between 
site i and % — 1 through the relation 

AUi = ln(fc/pi_a), (1) 



the average decrease of the potential per lattice spacing F = —hx{q/p) can be regarded 
as an average force which the particles are subjected to. Here and in the rest of the 
paper, the overbar stands for the average over p{p) and 7r(g). Without loss of generality, 
we assume F > 0, i.e. there is either an average bias to the right or the system is 
unbiased and we restrict ourselves to randomness distributions where the fraction of 
sites with pj < qj is finite. The natural control parameter /i of the one-particle problem, 
which we retain also for the many-particle process, is given by the positive root of the 
equation (for F > 0)[4J: 
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t =1- (2) 

It is monotonously increasing with F and zero in the unbiased case (F — 0). In the 

numerical simulations, we used a bimodal distribution for the hop rates, where = r 
holds for all i and the distribution of forward hop rates Pi is given by 

P (p) = c5(p-r) + (l-c)6(p-l), (3) 

with the parameters < r < 1 and < c < 1/2. For this distribution, the average bias 
is zero if c = 1/2 and the control parameter reads as 

In(Vc-l) (A] 

» = Hi/r) ■ (4) 



3. Steady state of a closed string 

We start our investigations by analysing the stationary properties of a finite closed 
string. In general, we are interested in the scaling behaviour of various quantities in the 
large L limit, when the global density of particles N/L is kept constant. 

With the purpose of studying the model in the presence of an external bias (// > 0), 
we invoke the random trap model of the random walk in a biased random-force type 
environment. This picture is based on the observation that the walker spends long times 
in certain localised regions (trapping regions) and in between, it performs a more or less 
directed motion. The process is thus approximated by a directed walk between traps 
characterised by effective trapping (or waiting) times, the distribution of which is broad 
[52] . This simplified model proved to describe the large-scale properties of the system 
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Figure 1. Illustration of the potential landscape denned by eq. |T]). The trapping 
regions are indicated by thick lines and the arrows show the heights of the barriers. 

correctly [I]. In the potential landscape of the infinite system, defined by eq. ([1]), which 
is a descending random walk path with average slope — F (see Fig. [I]), a trapping region 
can be identified as a basin on the left-hand side of a local maximum at site n for which 
U n > Ui holds for all sites i > n. The ascending section of the basin from the minimum 
to the maximum will be termed barrier. The waiting time r in a trapping region is in the 
order of the inverse of the persistence probability of a walker starting at the minimum 
of the basin with an imaginary absorbing site on its left-hand side [19J. This quantity 
is of the form of a Kesten- variable [33], the distribution of which has an algebraic tail 
for large r: 

p(r) ~ r- 1 -", (5) 

where /i is given by eq. ([2]). The waiting time r is related to the potential difference 
between the maximum and minimum, i.e. the height of the barrier U as r ~ e u for large 
U. For the ASEP with site-wise disorder, which has been studied in the framework of 
this simplified model [19], a phase separation can be observed: At the largest barrier, 
which serves as a bottleneck, a front appears, which separates a macroscopic high- 
density domain from a macroscopic low-density one and, as a consequence of particle- 
hole symmetry, the front is located where the potential is half the height of the barrier 
[34] . For the process under study, the situation is similar: the current is controlled 
by the largest barrier present in the system, where a front develops, however, as the 
particle-hole symmetry is broken, half-filling is no longer valid. Therefore, we consider 
first an isolated barrier of height U and express the current in terms of the one-particle 
waiting time. 

Let us assume that the front is located at some potential bU (b < 1), which is 
measured from the bottom of the barrier (see Fig. [2]). On the left-hand side of the 
front the particles sit closely next to each other so that the particle density is nearly 
1 (high-density phase), whereas on the right-hand side, particles are as far from each 
other as possible and the density is close to 1/(1 + 1) (low-density phase). In case of 
such a configuration, two kinds of processes which result in the shift of the front may 
take place: Starting from the front, holes diffuse to the left through the high-density 
phase and particles diffuse from the front to the right through the low-density phase 
over the barrier. In fact, the event that a hole detaches from the front and goes away 
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Figure 2. Typical configuration of particles at an isolated barrier for I = 1. 



from it is very rare: the probability for this decreases exponentially with the potential 
difference. Therefore the diffusing holes are practically alone in the high-density phase 
far from the front and the typical time scale for a hole to successfully overcome the 
potential barrier bU and escape to the left is in the order of e bU . The current of holes 
in the high-density region is thus Jh ~ e~ hU . The mechanism of particle transport in 
the low-density phase is similar. Here, particles are immobile since the interparticle 
spacings are almost everywhere maximal (hi = I). If the particle at the front jumps to 
the right, a short interparticle spacing with hi = I — 1 arises in front of this particle 
and thereby makes it possible for the next particle on its right-hand side, to jump to 
the right. If the latter event occurs, it can be regarded as if the short spacing moved to 
the right. For similar reasons as for the holes in the high-density phase, a short spacing 
is typically by oneself in the low-density phase and moves like a random walker taking 
steps of length of 1 + 1 lattice spacings. Thus, in order to escape to the right, this walker 
has to overcome only a reduced potential barrier (I — b)U/(l + I). The corresponding 
time scale is e*- 1-6 ^^ 1 ) and the current of particles through the low-density phase is 
Jp ~ i^j-e^ 1-6 ^/^ 1 ) . In the steady state, Jh = J P must hold, otherwise the front would 
move with a finite velocity. This leads to e~ bU ~ ^-L- e -( 1 - 6 ) c/ /('+ 1 ) an d for large U, i.e. 
U >> l\nl, we obtain b = 1/(1 + 2). The particle current j at the barrier is thus given 
by 

in terms of the one-particle waiting time r. 
3. 1 . Current 

It is clear that the quantity j derived above is merely the transport capacity of the 
barrier, i.e. the maximal current which can flow through it and the actual current is 
controlled by the incoming current of particles ji n if ji n < j. In a disordered ring of size 
L, the stationary current Jl is determined by the smallest one among the capacities of 
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barriers present in the system: J L = min{jj}. Making use of eq. ([SD and eq. (JSl), we 
obtain that the distribution of capacities has a power-law tail p(j) ~ j- 1 +( l + 2 )t i f or sm all 
j. The distribution of the sample-dependent current Jl is thus given by the well-known 
Frechet distribution for large L |35j : 

p( J) = (J + ^JW- 1 ^, (7) 

i 

in terms of the scaling variable J = cJ^L ( i+2 )f , where the constant c is related to the 
pre-factor in the asymptotical form of p(j). Thus, the current scales with the system 
size L as 

J L ~ l~ttt^ } (8) 

and vanishes in the limit L — > oo. 

We have performed numerical simulations for finite rings of size L = 
128, 256, 512, 1024, 2048 and after waiting sufficiently long time such that the system has 
settled in a steady state we measured the current. This procedure was then repeated 
for 3 x 10 4 independent samples for each L. As can be seen in Fig. [3] and Fig. HJ the 
Frechet distribution fits satisfactorily to the distribution of the current, and keeping in 
mind that the size of the largest trapping region is only O(lnL), the deviations can be 
assigned to corrections to scaling, which may be still considerable for the numerically 
available systems sizes. 




Figure 3. a) Distribution of the logarithm of the current measured in numerical 
simulations for different system sizes in the model with 1 = 2. b) Scaling plot of 
the distributions. The number of particles was N = L/2 and the binary randomness 
defined in eq. was used with c = 0.3 and r = 0.5, where the control parameter is 
fi 1.222. The solid curve is the Frechet distribution given in eq. ([7|). 



Reducing the control parameter, the exponent governing the finite-size scaling of 
the current is increasing and finally diverges as the unbiased situation is approached, 
i.e. when // — > 0. Strictly for /x = 0, the random trap approximation breaks down since 
the size of the largest trapping region is O(L). Since the height of the largest barrier is 
0(\/L), we expect that the magnitude of the current scales typically as — In | J L \ ~ a/Z. 
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Figure 4. a) Distribution of the logarithm of the current measured in numerical 
simulations for different system sizes in the model with 1 — 2. b) Scaling plot of 
the distributions. The number of particles was N = L/2 and the binary randomness 
defined in cq. ([3]) was used with c = 0.3 and r — 0.2, where the control parameter is 
fi ps 0.526. The solid curve is the Frechet distribution given in eq. ([7]). 



As the relaxation time is exponentially large, the dynamics in the unbiased case will be 
tested numerically in the context of the diffusion of a semi-infinite string in Sec. 

3.2. Density profile 

The steady-state average of the local occupation number Ui, which is zero (one) for 
empty (occupied) sites, is plotted against the site index i in Fig. [5] for a given sample. 
As can be seen, the steady state is segregated: the profile consists of a low-density 
phase, which extends from site ~ 150 to site ~ 320, and a high-density phase in the 
remaining part of the ring. In the low-density phase, the string is stretched out (i.e. 
hi = I) almost everywhere and this configuration is more or less fixed by the rugged 
landscape, which results in strong spatial variations in the profile with a period I + 1. 
After filtering out this noise (Fig. [51 lower panel), peaks can be observed in the profile, 
where the local density significantly exceeds the background density 1/(1 + 1). Similarly, 
in the high-density phase gaps with low local density develop. In the following, we shall 
estimate the number of peaks and gaps and will see that they both grow algebraically 
with L, however, with different exponents. 

The peaks and gaps are accumulations of particles and holes, respectively, at 
barriers with sufficiently large one-particle waiting times. First, let us consider a 
barrier in the low-density phase with a trapping time r and assume that the string 
is completely stretched out here, i.e. the interparticle spacings are maximal. Recalling 
the considerations of the previous section, the capacity of the barrier in case of this 
configuration is in the order of j^t^ 1 ^ 1 ^. Apparently, this configuration remains stable 
only if the incoming flow of particles Jl, which is determined by the largest barrier, does 
not exceed the capacity of the barrier concerning this configuration. In the opposite case, 
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Figure 5. Upper panel: Steady-state density profile of the model with I = 1 obtained 
by numerical simulation in a sample of size L = 400. The number of particles is 
N = 300 and the rates were drawn from the bimodal distribution with c = 1/3 and 
r = 1/4. Lower panel: Smoothed profile, obtained by averaging over pairs of sites: 
Pi = ((wi) + (^i+i))/2- The thick line is the potential landscape defined by eq. |T]) and 
the arrow indicates the starting point of the largest barrier. 



particles gradually accumulate here and form a high- density cluster. This process goes 
on until the front rises up to a certain potential level, so that the capacity concerning 
this new configuration is identical to the current Jl- (In this case, the magnitude 
of the potential at the front measured from the top of the barrier is jt^U, where 
U ~ (^+2) ln(l/ Jl) is the height of the largest barrier.) We thus expect an accumulation 
of particles at those trapping regions where j^jT^ 1 ^ 1 ^ < Jl holds, or, r > J L for 
large L. The number of barriers of this property is Ni ~ L J^-a+i) pMdr ~ LJl , 
where we made use of eq. (E]). Using eq. (jHJ), we obtain for the number of peaks in the 
low-density phase: 

Ni ~ LTT2. (9) 

In the high-density phase, the transport is realized by the diffusion of vacancies, 
and an accumulation, i.e. formation of clusters of I holes separated by a particle occurs 
at barriers for which r > J^ 1 holds. By a calculation analogous to that carried out for 
the peaks, we obtain that the number of gaps in the high-density phase is 

N h ~lM. (10) 

As can be seen, the exponents in eq. and eq. fllOp depend only on the parameter I and 
are independent of \l. In order to check these findings, we have performed numerical 



Anomalous transport 10 

simulations with a global density N/L = 2 (i+i) ' so b°th the length of the low- 
density phase and of the high-density phase was 0(L/2). It was measured in both 
phases how many times the smoothed steady-state profile intersects the horizontal line 
at the density 2 a±i) • ^ ne resur ts shown in Fig. [6] are in satisfactory agreement with the 
theoretical predictions. 




Figure 6. Size dependence of the average number of intersections (defined in the text) 
obtained by numerical simulation for (a) I — 1 and (b) 1 = 2. Bimodal randomness 
was used with c = 1/3 and r = 1/4 and the averaging was performed over 10 3 samples. 
The slope of the solid and dashed lines is 773 an( i j+k 1 respectively. 



3.3. Active particles, active holes 

In a given time in the steady state, some of the particles reside in the high-density phase 
and some of them reside in one of the particle clusters (peaks) in the low-density phase. 
Almost all of these particles are blocked by the hard-core exclusion. Other particles are 
located in the low-density phase and are blocked since the string is stretched out. All 
these particles may be termed inactive as they do not contribute to the current. On 
the other hand, there is a vanishing fraction of particles located between peaks in the 
low-density phase, which are not blocked by adjacent particles and which are responsible 
for the transport. In fact, the transport is realized by the diffusion of short interparticle 
spacings, which we shall call active particles. Analogously, in the high-density phase, the 
current is ascribed to vacancies which diffuse freely between gaps, the so called active 
holes. In the following, we shall estimate the number of active particles and active holes 
present in the system. 

In the high-density phase, the typical distance between adjacent gaps is ^ ~ 
L/Nh ~ L 1 ^ 5 . We will see a posteriori that the concentration of holes in these domains 
is vanishing, so that they can be regarded as independent walkers. It is known that the 
motion of a single random walker is controlled by the parameter \i If \i > 1, the 
velocity of active holes is constant, so they pass through the inter-gap region in a time 
t ~ C,h- The number of active holes N% present in this domain is related to the current 
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as N% ~ Ji,t. Using eq. ([HD, we obtain for the number of active holes between adjacent 

l-l/M 

gaps: N% ~ L l + 2 . If /i < 1, the diffusion of holes is anomalous, i.e. the travelling time 

is t ~ and the number of active holes between adjacent gaps is independent of L: 
N% ~ Jl^H^ ~ C(l). The concentration of active holes iV£/£/, is thus indeed vanishing 
for any /x. 

In the low-density phase, the typical distance between neighbouring peaks is 

a+i 

£j ~ L/A^ ~ L ; + 2 . The active particles (the short spacings) take steps of length I + 1, 
therefore their effective waiting time at a barrier is t 1 '" +1 ) . In the asymptotical form of 
the distribution of this quantity, (I + l)/i appears in place of /x, thus, we conclude that 
the motion of active particles is governed by the effective control parameter (I + If 
ji > the velocity of active particles is finite and we obtain for their number between 

two adjacent peaks: AT* ~ J L £i ~ L^~VT^. If ^ < i we have i ~ and the 

i 

number of active particles is independent of L: ~ Jl£i ~ 0(1). As can be seen 
also the concentration of active particles Ng/fy is vanishing. 

We have obtained that, as opposed to the number of peaks and gaps, the numbers 
of active holes and particles depend also on /i and the behaviour of the former changes 
at n — 1 while the behaviour of the latter changes at \i = j—. 

4. General formulation of steady-state results 

The results obtained in the previous section are easy to formulate generally for 
exclusion processes with random-force type disorder, where the active particles and 
holes responsible for the transport in the low-density and high-density phase (which 
are not necessarily true particles or holes, as we have seen above) overcome barriers 
according to the Arrhenius type activated dynamics characterised by a waiting time 
f. ^ e fiU _ Here, % = 1, 2 and the index 1 (2) refers to active holes (particles) in the 
high-density (low-density) phase, and the factors are characteristic of the particular 
system. This class of models includes the ASEP, the closed string, the ASEP with large 
particles [19| , the model with finite open strings to be discussed in Sec. O and possible 
combinations of these models. 

For models of this class, all the exponents characterising the stationary behaviour 
of quantities studied in the previous section can be expressed in terms of two exponents: 
lii = fi/ fi and /i 2 = /V/2, which appear in the distribution of effective waiting times fj 
of active particles and holes: Pi(fj) ~ Tj~ , i = 1,2. Generalising the argumentations 
presented in the previous section, one can easily show that the steady state current 

follows a Frechet distribution and vanishes according to 

1 

Jl ~ L . (11) 

The number of gaps in the high-density phase (i = 1) and the number of peaks in the 
low-density phase (i — 2) scales with the system size as 

1 M» 
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The typical number of active particles and holes in an inter-barrier domain is 




(13) 



For the closed string, the basic exponents are fii = \i and fi 2 = (I + l)fjt; for the ASEP, 
Hx = \i% = \i\ and for the exclusion process with particles of size d, which was studied 
in Ref. [19], \i\ = d\i and \ii = fi. 

5. Diffusion of a semi-infinite string 

So far, we have studied the steady state properties of the model defined on a ring. Now, 
we examine how a semi-infinite string moves when it enters an initially empty semi- 
infinite lattice (see Fig. [7]). To be concrete, the string enters at the first lattice site, 
which means formally that particles are created here with rate 1 and the particle closest 
to the entrance site is not allowed to jump to the right if it resides at site I + 1. We are 
interested in the time-dependence of the position (x) of the first particle (the head) of 
the string. Here, the angular brackets denote average over stochastic histories. Let us 



consider first an isolated barrier of height U and examine the way the head of the string 
goes over it. 

First, the front of the high-density cluster of accumulated particles is advancing 
(Fig. [5JA.) until the head reaches the top of the barrier (Fig. [SfB). This process is 
realized by diffusion of vacancies through the high-density cluster and its time scale is 
dominated by the escape time of vacancies when the head is in the vicinity of the top, 
which is 0(e u ). After the top was reached the head falls over and a stretched segment 
develops while the front separating it from the high-density segment moves backwards 
(Fig. [Hp). This process, which is realized by the escape of active particles through the 
stretched segment, goes on until the stable level of the front determined by the incoming 
current is reached (Fig. [Hp). Denoting the potential difference between this level and 
the maximum by Uo, the time scale of this process is 0(e u °^ l+1 ^). While for I > 1, the 
motion of the head is more or less directed after it has fallen over the maximum, for I = 1 
a special phenomenon can be observed: The head performs an unbiased random walk 
and may even return to the top (see Fig. [9]). As can be seen from the configurations 
shown in the figure, the stretched segment is located symmetrically on the two sides 
of the barrier hence the active particles go over the barrier in both directions with the 
same rate. Thus, the motion of the head can be described by a symmetric random walk 
with hop rates p(U') = q(U') ~ e~ u '^ l+1 \ which depend on the magnitude of potential 



co 




Figure 7. Advance of a semi-infinite string. 
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Figure 8. The mechanism of overcoming a barrier. 

U' (measured from the maximum) at the location of the head. The characteristic time 
in which the head first reaches a site of potential U' is in the order of e u '^ l+1 \ see e.g. 
Ref. [32] • However, during the swinging of the head, vacancies may also escape to the 
left through the high-density segment. These processes result in the shift of the front 
and thus induce a bias to the right for the motion of the head. The swinging therefore 
goes on practically until the front first reaches its stable level and afterwards the motion 
of the head becomes directed. 




Figure 9. The "swinging" of head of the string in the case 1 = 1. 



We see that the time-determining one among the above steps is the climb of the 
head to the top of the barrier, the characteristic time of which is in the order of the one- 
particle waiting time e . It is, however, valid only if the incoming current of particles ji n 
exceeds r _1 . Otherwise the time scale of overcoming the barrier is obviously 0(£/ji n ), 
where £ is the size of the trapping region. 

Now, we return to the investigation of the disordered model, where the potential 
landscape contains barriers of various heights. In what follows, we number the barriers 
starting from the first site of the lattice and measure distances in terms of the barrier 
index n. This can be done since the average distance between the starting points of 
consecutive barriers is finite. It is clear that ji(l + 1) is an upper bound for the velocity 
of the head, where ji is the capacity of the largest barrier the head has already left 
behind. This upper bound obviously changes at such a barrier which is larger than any 
of the barriers on its left-hand side. The mth barrier with this property, the so-called 
mth limiting barrier has a typical distance from the origin n m ~ 0(2 m ) and a typical 
trapping time r m ~ m^. Let us consider the mth limiting barrier, assume that m ^> 1 
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and estimate the time t that elapses until, starting from here, the head reaches the next 
limiting barrier, which is located a distance of 0(n m ) away. As aforesaid, the current of 
the limiting barrier provides a lower bound Tj for the travelling time: t > Tj ~ -. Tf^jr - 
On the other hand, the waiting time of the head at a barrier is at least in the order of the 
one-particle waiting time, therefore another lower bound is given by the single-particle 
travelling time: t > T x ~ J2i T i- Since the delay of the head at trapping regions and 
the inflow of particles at the mth limiting barrier take place simultaneously, we may 
write Tj,Ti < t < Tj + T\ and conclude that the travelling time is composed of two 
contributions: 

*~ 0(7)) +0(71). (14) 

i-i i 

The first term in eq. (1141) scales with the distance n m as Tj ~ n m {l+2)l1 ■ The second 



term is proportional to n m for fj, > 1, while it is 7\ ~ n^fj 1 for /i < 1. Comparing the 
two contributions, we obtain that, for /x > /x* = the travelling time is dominated 
by Tj, whereas for /i < /x* it is dominated by T\ for large n m . Therefore, inverting these 
relations, we obtain that the head of the string advances for /x > /x* asymptotically as 

(1+2) p 

(x) ~ t( ; + 2 ^+ 1 , (15) 

while if /x < /x*, we have 

(z) ~ t" (16) 

for large t. The diffusion of the head is thus anomalous and the diffusion exponent varies 

continuously with \i. At /i = \i\ a change occurs in the dependence of the diffusion 

exponent on /i and below this value, the motion of the head follows the diffusion law 

characteristic of the one-particle problem. 

The motion of the head between two limiting barriers is much more complex 

on a microscopic length scale. At small trapping regions with a waiting time Tj < 

3m ~ T m 1 i the head is delayed until the trapping region is filled up with particles. 

The release time of the head at large barriers with Tj > j" 1 , the number of which is 

n m /— i r _M- " 1 dr ~ n m r m''" i+2 ' ~ n^' 1 '^ l+2 \ is 0(rj). During this time, particles pile 

i_ 

up and a queue of length ^ ~ r m 1+2 Tj forms behind such a barrier. After the head was 
released, its velocity is temporarily determined by the capacity of the barrier instead 
of that of the true limiting barrier until the (the excess part of the) queue dissolves, 
which takes a time d/ji ~ (t^-) 1 ^^. For the length of the queue behind the m + 1st 

limiting barrier, we obtain £„ m+1 ~ r^ +1 ^^ +2 ^ ~ rim M - Thus for jj, > /x*, the length of 
the queue piling up behind a limiting barrier until the head is released is typically only a 
vanishing fraction of the distance from the preceding limiting barrier, while for /i < /i*, 
the domain behind the head is almost completely filled up with particles. In fact, in 
the latter case, also queues forming at non-limiting but large barriers may extend to the 
preceding limiting barrier. The queue temporarily blocks the inflow of particles here, 
which explains the increase of travelling time of the head compared to that dictated 
solely by the current of the limiting barrier. We see again from the dynamics of the 
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queues forming at trapping regions that the properties of the system change at the value 
//* of the control parameter. 

We have performed numerical simulations and measured the time-dependence of the 
average displacement of the head for different values of the control parameter. Results 
are shown in Fig. [10J As can be seen, for fi > //*, the diffusion exponent is in good 
agreement with the predictions. For and for short times, the diffusion seems 

to follow the law found for fi > fx* (eq. (ITS]) ) and for longer times a crossover can be 
observed to the true asymptotic behaviour (eq. ( |T6i) ). This means that for short times, 
the second term in eq. (JHj) is suppressed by the first one, but ultimately, it dominates 
the travelling time since it grows with a greater exponent. 




Figure 10. Results of numerical simulations. Left panel: the logarithm of the average 
position of the head plotted against the logarithm of time in the model with I — 1 for 
different values of the control parameter. Binary randomness was used with parameters 
c = 0.2, r = 0.5; c = 0.2, r = 1/8; c = 1/3, r = 1/8; c = 1/3, r = 1/16, where the 
control parameter takes the values /j, = 2, 2/3, 1/3 and 1/4, respectively. The average 
was performed over a few hundred independent samples and 3 — 20 runs for each 
sample. The slope of the solid lines is > w ^i\e that of dashed lines is p. Right 

panel: Square root of the average position of the head plotted against the logarithm 
of time in the model with I = 1 in the unbiased case /i = 0. The parameters c = 0.5 
and r = 0.5, 0.2 were used and the data were averaged over 200 independent samples 
and 20 runs per sample. 



In the limit /i — > 0, the diffusion exponent tends to zero and we expect that for 
/i = 0, the displacement (x)(t) grows slower than any power of t. The unbiased case 
is out of the scope of the above theory, nevertheless, we can estimate the leading order 
behaviour of (x)(t) by taking into account that the sojourn time t of a random walker in 
a domain of size £ scales as hit ~ Thus, we expect that the average displacement 
grows as (x)(t) ~ C(lnt) 2 , where C is a non-universal constant. Results of numerical 
simulations are in accordance with this relation (see Fig. [TU]) . 
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In the final part of this work, we shall study finite open strings composed of m particles, 
which we call shortly m-strings. 

6.1. Steady state 

First, we are interested in the steady-state properties of systems defined on a ring, in 
which the density of m-strings is finite. This model belongs to the model class discussed 
in Sec. HI for which the steady-state results presented there are at our disposal once the 
basic exponents \L\ and [i^ are known. 

In the high-density phase, the active holes which are responsible for the current in 
the inter-barrier regions are simple holes, consequently, we have \x\ = fi. In the low- 
density phase, the active particles are m-strings moving in the inter-barrier domains. 
Thus we need the diffusion exponent of an m-string characterised by the parameters 
m and I. This can be obtained by constructing the network of possible transitions of 
such systems. Let us consider the simplest case, m = 2, a "dimer" and denote the state 
when the left-hand side particle of the dimer is located at site % and the right-hand 
side particle at site i + 1, i + 2, i + 3, . . . by A iy B i} C i: . . ., respectively. The network of 
transitions for I — 3 is depicted in Fig. [IT] and the corresponding networks for I = 2 
(I = 1) are obtained by deleting the nodes Di (Di and Cj) and the links connected with 
them. The conclusions drawn through this example can be easily extended to general 




Figure 11. Network of transitions for an m-string with m = 2 and I = 3. An arrow 
with index i symbolises an allowed transition in the direction of the arrow with rate 
Pi and in the opposite direction with rate qi+i- 

m and I. As can be seen, the diffusion of an m-string can be regarded as a random walk 
on a disordered quasi-one-dimensional lattice. The walker has several alternative paths 
to get from one site to another one, however, the potential difference defined by eq. (00) 
is the same along all paths in an infinite system. Thus, a single-valued potential can 
be introduced here, just as for a one-dimensional chain. Therefore, when getting from 
one site to a remote one, the trapping regions cannot be "walked round" and in the 
point of view of large-scale properties, the existence of alternative paths is irrelevant. 
Moreover, one can easily make sure that following a path in the network between two 
distant points, each type of link is contained m-times in the path. See, for example, the 
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path Ai — > -Bj — > — > -Bj+i — > . . . in the figure, in which the links of the original 
lattice are contained doubly. Consequently, the effective height of a barrier felt by an 
m-string is m times greater than the true height. This can also be seen in a more 
heuristic way by taking into account that one step of an m-string amounts to m steps 
of its constituent particles. Thus, we conclude that ji — m and the exponent governing 
the diffusion of an m-string is I- independent and is given by fi% = fijm. The results in 
eq. (ITT1) - ([TB1) are thus valid for this model with \x\ = fi and fi 2 = fi/m. For example, for 
the finite-size scaling of the current we obtain: 

J L ~ L ^ m + 1 . (17) 



6.2. Invasion 

In the following, a problem similar to that discussed in Sec. |5] will be examined, i.e. 
we consider a semi-infinite lattice, which is initially empty, and at the first lattice site 
m-strings enter ("invade") the lattice with a finite entrance rate and we are interested 
in the advance of the m-string first entered the system, as well as the evolution of 
the number of particles on the lattice. Although, the steady-state behaviour could be 
handled together with the problem of a closed string within a common formalism, the 
problem of invasion of m-strings must be treated separately from the diffusion of a 
semi-infinite string and it has much in common with the invasion of uncoupled particles 
(m = 1), which has been studied in Ref. [T9] . 

From now on, we shall call the m-strings particles. First, it is clear that the 
first particle advances with a constant velocity if // 2 > 1 {fi > m). Let us assume 
that < fi2 < 1 and consider the kth limiting barrier far from the entrance site, i.e. 
k 1, and estimate the travelling time of the first particle between the kth and the 
k + 1st limiting barrier. In this domain, the first particle is caught up by the following 
particles only at sufficiently large barriers, for which f 2 > j^ 1 , where jk ~ n k 1// ^ 1+M2 ^ 
is the current supplied by the kth limiting barrier. The number of such barriers is 
n B ~ n^ 1 ^ 14 " . Such a trapping region is filled up with particles in a period of O^ 1 ) 
to a level at which the escape time of the first particle is in the same order of magnitude, 
0{j k l ). The waiting time of the first particle at such a barrier is thus 0(r4 + ). 
The travelling time of the first particle between two such adjacent barriers, which are 
located a distance £ ~ n k /n B ~ n ^/(^+^) a p ar t i s i n the same order of magnitude: 
to ~ £o^ 2 ~ nl/^^^. Thus the total travelling time between the kth and the k + 1st 
limiting barrier is t ~ nBn]/^ 1+ ^ ~ n ^+^)/^+^) _ Inverting this relation, we obtain 
that the advance of the first particle follows the asymptotic law: 

l+l/m 

(x) ~ t 1 ^^, fi < m. (18) 

Nevertheless, the average number of particles N(t) can be shown to grow slower. 
The number of particles present in the system is approximately equal to the length of 
the high-density segment, which extends from the entrance site to an advancing front 
|19j . When the front is located at a limiting barrier in a distance iV from the entrance, 
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the inflow of particles is controlled by the capacity of this barrier. For the growth rate 
of N, we may thus write ^ ~ jV-VOi+z^ w hi cri yields 

1+1/ m 

N(t) ~ (19) 

So, the total number of particles in the system grows with a smaller exponent than the 
displacement of the first particle. Such a "dispersion" is obviously not possible in the 
case of a semi-infinite string since particles are coupled. Setting m = 1 in the above 
expressions we recover known results for the ASEP. 

In the unbiased case (/i = 0), the advance of the leading particle is expected to be 
ultra-slow and to follow the logarithmic scaling law characteristic of Sinai diffusion. 

7. Discussion 

As it was outlined in Ref. [19] . the holes in the closed string with I = 1 and with hop 
rates Pi,qi can be regarded as large particles of size d = 2 with hop rates pi = 
Qi — Vi-x- Although, for I > 1, a strict mapping between the two models does not exist, 
the phenomenological results obtained in this work show that the basic exponents fi2 
of the closed string with parameter / are related to the basic exponents p>i,p>2 °f an 
exclusion process with particles of size d — I + 1 in the following way: H\ — fa, 1^2 — fa- 

Comparing the exponent governing the decay of the steady-state current of a closed 
string to that of the ASEP, we see that the former is smaller for any /i > 0, i.e. 
the coupling between particles facilitates the transport in this set-up. Moreover, the 
exponent decreases monotonously with I and finally tends to zero in the limit I — > 
When not all particles are coupled with each other but they form finite m-strings, the 
coupling has the opposite effect: With increasing m, the transport slows down more and 
more compared to the ASEP both in the steady state and in the case of the invasion. 

Contrary to the problem of invasion of uncoupled particles, in the case of a semi- 
infinite string, both the displacement of the head and the number of particles in the 
system grows obviously with the same power of time. This exponent is larger than that 
describing the increase of the number of particles in the ASEP if fi > 1/2, meaning 
that the coupling is favourable in this case in the point of view of the bulk of particles. 
However, if /x < 1/2, the inflow of particles is slowed down by the coupling. As far as 
the first particle is concerned, it diffuses faster in the ASEP if /i > /i*, otherwise the 
diffusion exponents are equal in the two models. The influence of the parameter I on the 
transport shows a tendency similar to that in the closed string: The diffusion exponent 
increases monotonously with / (except of the regime //<//*, where it is /-independent) 
and tends to 1 if I — > oo. 

The diffusion of the semi-infinite string can be speeded up by "pulling" the head of 
the string, that means when the first particle is not allowed to hop to the left. In this 
case, one can show that the time scale of overcoming a barrier is reduced from 0(t) 
to 0{t x '( 1+2 >) } which results in that the diffusion law in eq. ( fl5|) is valid in the entire 

| We remind the reader that our results are valid in scaling regime InL ^> I. 
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biased phase /i > 0. Comparing this model to the invasion of uncoupled particles, the 
diffusion of the pulled string is faster than that of the first particle in case of the ASEP 
only if /i < j^, while the inflow of particles is faster for the pulled string for any /i. 

Relaxing the hard-core exclusion condition in the ASEP, i.e. prescribing for the 
positions of consecutive particles the weaker condition Xi < Xi + i instead of Xi < Xj+i, 
we arrive at a zero-range process. Doing so for the model with coupled particles but 
keeping the attractive interaction x i+ i — Xi < 1 + 1, we obtain a model which is no longer 
a ZRP. Due to the absence of hard-core repulsion, the high-density phase in the steady 
state of the closed string shrinks to the minimum in the largest trapping region, where a 
condensate forms and the low-density phase extends to the whole system. The current 
is thus expected to vanish as Jl ~ L" 1 ^ 2 ~ L~ < ;+1 )^ . 

Finally, we mention that, as a possible generalisation of the model with limited 
inter-particle spacings, one could apply elastic pair-interactions by the help of which 
the crossover between the model with uncoupled particles (ASEP) and that studied in 
the present work could be investigated. 
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